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Abstract: We compute one loop neutrino masses in the R-parity violating Minimal 
Supersymmetric Model, including the bilinear ^masses in the mass insertion approx- 
imation. To the order we calculate, our results are independent of the Higgs-lepton 
basis choice. We have a variety of perturbative parameters — gauge, yukawa and 
trilinear couplings, and R p violating masses. Their relative magnitudes determine 
which diagrams are relevant for neutrino mass calculations. We find new loop dia- 
grams which can be relevant and have frequently been neglected in the past. For the 
Grossman-Haber neutral loop contribution to the neutrino mass matrix we obtain 
explicit analytic results. 
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1. Introduction 

An interesting challenge in particle physics is to construct a consistent model which 
can accomodate neutrino oscillations in order to explain atmospheric and solar neu- 
trino data [1]. Considering the results of SuperKamiokande and all solar neutrinos 
experiments, at least two different values of Am 2 and sin 2 29 are necessary to account 
for the data. Neutrino masses cannot be generated in the Standard Model (SM), so 
it is necessary to consider extensions of the SM [2]. An interesting possibility for 
generating neutrino masses are models in which lepton number violation occurs. In 
these models Majorana neutrino masses are generated by interactions which violate 
lepton number in two units: AL = 2. 

Among the possible models we consider the Minimal Supersymmetric Stan- 
dard Model without imposing _R-parity. The quantum number R is defined as 
R p = (_i) L + 3B + 25 ' ) where L, B, S are the lepton and baryon number and the spin of 
the particle, respectively [3]. Imposing i?-parity in a supersymmetric model has the 
advantage of suppressing simultaneously the presence of baryon and lepton number 
violating interactions, which are strictly constrained from proton decay experiments 
and other low-energy processes [4]. However, one can impose a less stringent con- 
straint in which all interactions conserve baryon number. Thus only lepton number 
violating terms are present in the Lagrangian. These terms allow left-handed neu- 
trinos to obtain a Majorana mass, at tree level through mixing with the neutralinos, 
and through loop diagrams that violate lepton number (in two units) [5, 6, 7]. 

In the SM, the Higgs and leptons have the same gauge quantum numbers. How- 
ever, they cannot mix because the Higgs is a boson and the leptons are fermions. 
In a supersymmetric model this distinction is removed, so the down-type Higgs and 
sleptons can be assembled in a vector Lj = (H d , Lj) with J : 4..1. With this notation, 
the superpotential for the supersymmetric SM with R p violation can be written as 

W = fi J H u Lj + \ JKe LjL K E c e + \' J ^LjQ p D c q + h p t q H u Q p U c q (1.1) 

The R p violating and conserving coupling constants have been assembled into vectors 
and matrices in Lj space: we call the usual \x parameter /i 4 , and identify the usual 
6i = Hi, \h?^ = \ 4jk , and h^ 9 = \' 4pq . Lower case roman indices i,j,k and p,q are 
lepton and quark generation indices. We frequently suppress the capitalised indices, 
writing (I = (/i 4 , /i 3 , fj, 2 , fJ>i)- 

We also include possible R p violating couplings among the soft SUSY breaking 
parameters, which can be written as 

= ^HlH u + \L J \m\\ JK L K + B J H u Lj 

+A u v s H u Q p U c s + A j p s LjQ p D c s + A JKl LjL K Ef + h.c. . (1.2) 
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Note that we have absorbed the superpotential parameters into the A and B terms; 
e.g. we write B 4 H u H < i not B 4 fj J 4 H u Hd 1 . We abusively use capitals for superfields 
(as in (1.1)) and for their scalar components. 

We have put the Higgs into a vector with the sleptons, and combined the R p 
-violating with the R p conserving couplings, because the lepton number violation can 
be moved around the Lagrangian by judiciously choosing which linear combination 
of hypercharge = -1 doublets to identify as the Higgs/higgsino, with the remaining 
doublets being sleptons/leptons. However, this freedom to redefine what violates 
L is deceptive, because phenomenologically we know that the leptons are the mass 
eigenstate e, fi and r , so we know what lepton number violation is. Lepton number 
is defined in the charged lepton mass eigenstate basis — the freedom to choose which 
direction is the Higgs in the Lagrangian just means that there is not a unique inter- 
action eigenstate basis. There are two possible approaches to this fictitious freedom; 
choose to work in a Lagrangian basis that corresponds to the mass eigenstate basis 
of the leptons, or construct combinations of coupling constants that are indepen- 
dent of the basis choice to parametrise the R p violation in the Lagrangian [6, 8, 9]. 
These invariant measures of R p violation in the Lagrangian are analogous to Jarlskog 
invariants which parametrise CP violation. 

The standard option is to work in a basis that corresponds approximately to the 
mass eigenstate basis of the leptons. For instance, if one chooses the Higgs direction 
in Lj space to be parallel to fij, then the additional bilinears in the superpoten- 
tial Hi will be zero. In this basis, the sneutrino vevs are constrained to be small 
by the neutrino masses, so this is approximately the lepton mass eigenstate basis. 
Lepton number violation among the fermion tree-level masses in this basis is small 
by construction, so it makes sense to neglect the bilinear R p violation in setting 
constraints on the trilinears, as is commonly done (for a review, see .e.g. [4]. For a 
careful analysis including the bilinears, see [10]. Note that in many cases, the most 
stringent constraints on the additional R— parity violating parameters come from 
neutrino anomaly data [11].) 

The aim of the "basis-independent" approach is to construct combinations of cou- 
pling constants that are invariant under rotations in Lj space, in terms of which one 
can express physical observables. By judiciously combining coupling constants one 
can find "invariants" which are zero if R p is conserved, so these invariants parametrise 
Rp violation in a basis-independent way. For instance, consider the superpotential 
of equation (1.1) in the one generation limit, I : 4. .3. It appears to have two Rp 
violating interactions: {i 3 H u L and \'LQD C . It is well known that one of these can 



1 We do this because Bj is a vector — a one index object — in {Lj} space. From this perspective, 
giving it two indices can lead to confusion. 
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be rotated into the other by mixing H d and L [5]. If 



L' = -fJ^H.-^J^L , (1.3) 

then the Lagrangian expressed in terms of H' d and U contains no H U U term. One 
could instead dispose of the \'LQD C term. The coupling constant combination that 
is invariant under basis redefinitions in (H d , L) space, zero if R parity is conserved, 
and non-zero if it is not is ji^X' — h d fi 3 = (/i 4 , /x 3 ) A (h^X'). For a more detailed 
discussion of the approach followed in this paper see eg [8, 9, 12]. 

Basis-independent coupling constant combinations which parametrise the amount 
of Iflp can be constructed in various ways. Their advantage over the coupling con- 
stants is that they cannot be set to zero by a basis transformation. So we define the 
following invariants {5} which reduce to the coupling constant of the same indices 
in the basis where the sneutrinos do not have vevs (the — > is to what the 5 becomes 
in this basis). 

(1.4) 

\'PQ . \i . v 

4 = 4^-^ ("J 



\B\mf B. 
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We have chosen the Q and D c bases to make v ■ \' pg diagonal — that is, to diagonalise 
the down-type quark mass matrix. We assume that this is also the mass eigenstate 
basis for the squarks. Similarly, we choose the E c basis to make v ■ X t X k ■ v oc 5 lk . 
Combined with our definition of the lepton directions as Lj = v-X l /m^, (mf = \v- X 1 \) 
this means we diagonalise the charged lepton mass matrix induced by the Higgs 
vev. Here we are neglecting i?-parity violating bilinears that mix the charginos and 
charged leptons; these masses will be included as mass insertions in perturbation 
theory. 

The mass matrix of the neutralinos is composed by the 4 neutralino fields of the 
R-parity conservng MSSM and the 3 neutrino fields. In an interaction eigenstate 
basis they are 

Xint = (B, W 3 , h u , h d , v T , Vp, u e ) (1.8) 

where as previously discussed, in fi p models there is no unique interaction eigen- 
state choice for the basis in hd and lepton space. In this paper we define the basis 
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where the sneutrinos do not have vacuum expectation values to be the interaction 
eigenstate basis. We number the elements of x m reverse order, so xJnt = B, and 
) = {h d ,v T ,v ll ,v e ). We use this numbering so that (x 3 , X 2 , X 1 ) will 
be the mass eigenstate neutrinos. 

In the R p conserving MSSM, lepton number is conserved, and the first four 
fermions have majorana masses, while the three neutrinos are massless. The mass 
matrix can be diagonalised with a matrix : 



where Z is a 4 x 4 matrix and / is the 3x3 identity matrix. 

The purpose of this paper is to address and clarify the issue of different basis in 
ftp neutrino mass calculations. We construct "basis-independent" estimates of each 
loop diagram contributing to the neutrino mass matrix. These estimates are in an 
arbitrary Lagrangian basis. They lead us to include the bilinear ^masses in the mass 
insertion approximation. This introduces a new diagram, and new contributions from 
the usual diagrams. These new terms resolve the basis related puzzles that arise in 
ftp neutrino mass calculations when the bilinears are neglected in the loops. Clearly 
neutrino masses should not depend on the basis they are computed in; however, if 
the ftp masses are neglected in the loops, then what one is neglecting depends on 
the basis. We will return to this in the discussion at the end. Having verified the 
irrelevance of performing basis independent calculations, we present estimates and 
the calculation of the neutral loop in the basis where the sneutrinos do not have 
vacuum expectation values, because the results are more compact than the basis 
independent estimates. 

In section 2 we introduce our procedure for obtaining neutrino masses, and in 
section 3 we estimate the basis-independent contributions from all loops. In section 
4 we provide complete analytic expressions for the neutral loop, and compare with 
the exact result previously obtained in the literature [12]. The latter calculation is 
done in the MSSM mass eigenstate basis. The neutral graph is of interest as it 
has been analysed very little in the literature and furthermore it can provide strong 
constraints on Higgs physics in this model. In section 5 we discuss in more detail the 
issues related to the basis choice and conclude. 

2. Perturbation theory with fo p masses 

Neutrino masses are observationally known to be small, so we can compute them 
in perturbation theory. There are two ways to include the ft p bilinears in one-loop 
calculations of neutrino masses. One approach is to diagonalise the tree level mass 
matrices for neutral and charged scalars and fermions, and calculate the one-loop 
contributions to the neutrino mass matrix using a tree-level mass eigenstate basis 
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for the propagators in the loop [13]. This is usually done numerically. Note that for 
5^ 7^ 0, 5b 7^ 0, the charged and neutral Higgses [MSSM neutralinos] mix with the 
sleptons [leptons]. The alternative, which we follow here, is to propagate particles 
in the MSSM mass eigenstate basis, and include the ft p masses as mass insertions 
[14] in perturbation theory. For "sensible" basis choices that are close to the MSSM, 
the 1/lp bilinear masses and couplings are small, so including them as mass insertions 
should be adequate. 

The algorithm for computing neutrino masses, with Ifi p bilinears included in the 
mass insertion approximation, consists of the following steps: 

1. choose a basis where the R p violating parameters are small, eg the basis where 
the sneutrinos have no vevs, or where there are no Hih u £i terms. 

2. diagonalise the heavy neutralino 4x4 matrix. 

3. calculate the tree-level seesaw neutrino mass. 

4. calculate loop contributions to the 3x3 neutrino mass matrix, including the 
]/{ p bilinear masses as mass insertions. 

5. diagonalise the resulting 3x3 neutrino mass matrix. 

It is well known that in the ]fl p MSSM only one of the neutrinos acquires a 
mass at tree-level through the see-saw mechanism due to the mixing of the neutrinos 
and neutralinos. We will focus on the (finite) loop contributions to neutrinos which 
are massless at tree level. If the tree-level mass is non-zero, there are gauge loop 
corrections which must be renormalised. To avoid these, we for simplicity neglect 
loop corrections to the neutrino who is massive at tree level. We will address loops 
containing gauge bosons, and issues of renormalisation, in a subsequent publication. 
If the tree-level mass is smaller than the loop contribution, it can be set to zero 
(Sp = 0), and one- loop masses can be calculated for three neutrinos. 

We only consider one-loop diagrams contributing to the neutrino mass matrix- 
higher order gauge loops may be larger than one loop diagrams with Yukawas/trilinears 
at the vertices, but gauge couplings have no flavour indices so pure gauge loops will 
give a mass aligned with that of some lower order diagram. Two-loop diagrams with 
Higgs exchange rather than mass insertions on the internal lines should be suppressed 
by an additional factor of l/(167r 2 ). We neglect the loop corrections to the neutralino 
sector and to neutrino-neutralino mixing. These would contribute to the two loop 
neutrino mass matrix. To compute the one-loop neutrino mass matrix, we propagate 
tree level mass eigenstates in the loop with tree level interactions. 

In this paper we concentrate on step 4. We draw diagrams and make basis- 
independent estimates in an arbitrary Lagrangian basis, which we assume satisfies 
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step 1. The advantage of this is that the basis-independent coupling constant com- 
bination controlling the neutrino mass can be read off these diagrams. However, to 
get the correct dependence on MSSM masses and mixing angles, we should work 
in the MSSM mass eigenstate basis. We do this for the neutral loop in section 4, 
defining the Higgs to be the direction in Lj space that gets a vev. This basis choice 
in L j space is close to the charged lepton mass eigenstate basis because the mis- 
alignment between \x and v must be small, so /ij in this basis are small. There are 
also no lepton number violating D-term masses proportional to the sneutrino vev. 
Since we are perturbing in lepton number violating masses, this reduces the number 
of interactions we must consider. This is therefore a good basis from which to do 
perturbation theory in $ p parameters. We define our "MSSM mass eigenstate basis" 
(in an ^theory) by diagonalising the mass matrices with all the ^parameters set 
to zero, then reintroducing the ft p parameters in the Lj basis where the sneutrino 
have no vevs. 

We have many parameters in which we will be perturbing to compute loop neu- 
trino masses: Standard Model Yukawa couplings h and gauge couplings g, of very dif- 
ferent sizes, and bilinear and trilinear ft p couplings of unknown size. We parametrise 
the size of the ft p couplings in a basis-independent way via the parameters 5 intro- 
duced earlier. We would like to order them in magnitude, so that we consistently 
include all contributions greater than a certain size. 

It is clear that h <^ g, and that the Yukawas are known and of different sizes. 
We do not know the size of the various {5}; so we consider three possibilities 2 : 

• A: Assume bilinear ft p is small and can be neglected in the loops. This is the 
case if h e 5\, hdSy ^> <?<5b, gtf^. One gets the usual loop diagrams of figures la 
and lb. 

• B: Neglect <5 M , but include the soft bilinear ft p terms (Bi and m^) in perturba- 
tion theory. This means h e 5\,hd5\> ~ g5s, but g5^ <C h e 5\,hdSy. This gives 
the neutral Grossman-Haber loop of figure lc, and additional contributions to 
the usual loop of figure la. 

• C: Include all bilinear ft p terms in perturbation theory. There are additional 
contributions to diagrams la, lb, and lc, and a new diagram Id. 

Note that models with only bilinear ft p can fall into any of these three categories. 
A model where the misalignment between the vev and the Yukawa couplings is much 
greater than between the vev and jl would fall into case A. Case C corresponds to a 
model where the misalignment between the vev v and \x is similar to the misalignment 
between the Yukawa couplings and v. 

2 We have suppressed the family indices for simplicity. 
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3. Estimates 



The possible diagrams that contribute to the neutrino mass must contain two scalar- 
fermion-fermion couplings at the vertices of the loop, and two AL = 1 interactions. 
The loop can contain coloured and colour-singlet charged or neutral particles, and we 
can have either two gauge couplings, one gauge and one Yukawa/trilinear coupling, 
or two Yukawa/trilinear couplings at the vertices. 

The charged coloured loop is figure lb. This has Yukawa or trilinear couplings 
at the vertices, where a higgsino or neutrino interacts with a quark and a squark. 
Lepton number violation can be due to the trilinears, or to mass insertions mixing 
the neutrino with a neutralino on the external legs. 

In figure lc we have the contribution from the neutral loop. At the vertices we 
must have two gauge couplings. It is not possible to have vertices with combinations 
of gauge- Yukawa and Yukawa- Yukawa couplings for a neutral loop, because the vertex 
with the Yukawa coupling would involve an E c . This diagram has been previously 
considered in [7], and we will refer to it as the Grossman-Haber (GH) diagram. 
Lepton number violation, AL = 2, occurs due to mass insertions on the internal 
scalar propagator. In the < v >= basis, we cannot have a neutral loop with one 
unit of lepton number violation on each of the propagators because if we turn the z 
into a v, we would also need d^u — > Z M . This is not possible, because the Goldstone 
boson is in the direction of the vevs, so has no v component. 

There are in addition two charged loops, drawn in figures la and Id. Charged 
loops cannot have two gauge couplings at the incoming and outgoing vertices because 
AL = 2 on a charged line is forbidden by charge conservation. It is not possible to 
have AL = 1 on each line for the same reason as for the neutral gauge loop. One 
gauge vertex will produce a charged slepton and a wino/chargino; if w + — > t by 
a AL = 1 mass, then d^L + — > W 1 would be required on the bosonic line. The 
Goldstone boson has no L + component in the < v >= basis. For charged loops 
we can have two Yukawa/trilinear couplings at the vertices, which is the canonical 
trilinear diagram of figures la. 

It is important to note that charged loops can have one gauge and one Yukawa 
coupling at the vertices. We show this contribution in figure Id. To have one gauge 
and one Yukawa/trilinear coupling, we need (charged) gaugino-lepton mixing on the 
fermion line. This means we need 5^ ^ 0, so these diagrams only contribute when 
we make assumption C. We will discuss the supersymmetric partner of diagram Id, 
which could arise if mixes with d^E , in a subsequent publication. 

In the following subsections we make basis-independent estimates of the size of 
each loop, assuming that the Higgs vacuum expectation value, Higgs masses, slepton 
masses and chargino/neutralino masses are all of the same order ~ msusY- We 
will calculate the dependence of the neutral loop on the masses of the propagating 
particles in section 4. We also present estimates in the basis where the sneutrinos 
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do not have a vev, (u) = 0. A more detailed discussion of the diagrams involving 
charged loops will be given in a subsequent publication. 

3.1 5fj,, 5b ~ 0: the canonical loop 

Figures la and lb contain the usual diagrams considered in the literature contributing 
to one loop neutrino masses when we allow the lepton number violation to occur only 
at point II and VII. We neglect the Iji v bilinear masses in this section, so the only 
available sources of ^? p are the trilinear couplings at the vertices. Lines 1,2 of table 1 
list the basis-independent coupling constant combinations to which these diagrams 
are proportional and estimate the contribution to the neutrino mass matrix. Note 
that on both the scalar and the fermion lines the particle must flip between doublet 
and singlet via an interaction with the Higgs/slepton vev. We represent this as a 
mass insertion of the vev to assist in writing basis-independent estimates. 

To get the estimate in column 4 of table 1, we assume that we have some incident 
z/j, in some choice of basis. At vertex II, we have a trilinear coupling X lJk , where in 
a generic basis we allow J : 1..4 and k : 1..3. A slepton propagates to the top of 
the loop, where a trilinear mass insertion A RTk VR + X RTk fiRV u turns the Ef. into an 
Lt- We again generically allow T : 4..1 — this means we have not chosen the singlet 
charged lepton basis. L T propagates to the vertex VII. Now following the fermion 
line down, we propagate an £j to the bottom of the loop, flip it into an SU(2) singlet 
via the mass %A JM ", and propagate e c n to vertex VII, where we put a trilinear 

X Tjn 

Thus, the contribution to the neutrino mass matrix is 

\ ijk \Tjn 

my oc X JMn v M [{XA) RTk v R + v u( , R X RTk ], (3.1) 

107T m SUSY 

where all repeated indices are summed. In basis where sneutrinos do not have a vev, 
and assuming that (XA) RTk = A\ RTk 3 this is approximately 

Sinkdkn e e 
^ °\ °\ m n m k fon\ 

16tt 2 m SUSY 

where rrf n = \v • X n \ are the charged lepton masses generated by the Higgs vev, for 
n : 1..3. We can compare this expression to the more correct formula given by 

m .. „ {x ink X^ kn m e n mi(A + //tan/3)} f ^ k \ , (3.3) 



where L k = m LU /m E . 



/<*> = tH- < 3 - 4 ) 



3 This avoids flavour violation among the sleptons. Additional diagrams would appear when one 
allows flavour violation in the slepton sector [15]. 
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and \i= \\x\. 

Similarly, figure lb with A' couplings produces a neutrino mass matrix 

o \ 'isp \ 'jrt 

my oc -^-5 A' Ms ^ M ((AA)'^ + t^A'**). (3.5) 

i07r m SUSY 

Recall that we chose the quark basis to diagonalise the mass matrix [v ■ \'} pg . The 
exact result is 

my ~ 3X' ips X' jsp m d p m d s {A + p tan 0) J^ p \ , (3.6) 

J.D7T Tfl p>c 

v 

where Q p = rn 2 Q P /m 2 D c. 

3.2 <?„ ~ 0, 5 B ^ 0: the GH loop 

As mentioned above, there can be a loop contribution of order g 2 to the mass of 
the neutrinos that are massless at tree level from the diagram of figure lc. This has 
been discussed in detail be Grossman and Haber [7] (see also [16, 12]). The soft 
SUSY breaking terms B k and m\ k , and the sneutrino vev, mix the Higgses with the 
sleptons. If the soft masses are universal and B k = B/ik, diagram lc will be a loop 
correction to the mass of the neutrino that is massive at tree level. In figure 3 we 
show the supersymmetric partner of this diagram which gives a small contribution 
to the mass of the neutrino which is massive at tree level. However if B is misaligned 
with respect to both v and jl, the diagram of figure lc will give a mass also to the 
neutrinos that are massless at tree level. 

The diagram in fig.lc is in the "MSSM mass eigenstate basis" , where the sneutrino 
does not have a vev. The sneutrinos mix with the neutral CP-even Higgses h, H 
and the CP-odd Higgs A through B k and m\ k . In Section 4, we calculate this 
diagram in the mass insertion approximation, and show that we get the same answer 
as by expanding the exact result of [12]. We also show the contributions to this 
diagram in an arbitrary basis in figure 2, because this facilitates basis-independent 
estimates. There is an additional diagram (not drawn) combining those of figure 2, 
with one unit of lepton number violation from B and one from fh 2 . In an arbitrary 
interaction eigenstate basis, the sneutrino could have a vev, which induces lepton 
number violating masses through the D-terms. This would correspond to figure 2 
with diagonal fh 2 and L violation from the sneutrino vevs. 

Figure 2 and the diagrams with one unit of lepton number violation from B and 
one from fh 2 contribute to the neutrino mass matrix by 

rriij oc g 2 (v u Bi + m 2 K v K )(v u B j + v M m 2 Mj ). (3.7) 

We can write f jfrn^]^ in terms of Bi using the minimisation conditions for the 
Higgs potential [12]: 

v u Bj + {5 JK m 2 z cos 2(3 + [m L ] 2 JK )v K = 0, (3.8) 
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where at tree level [mt]|j = m 2 JK + /ij/ir, and B and m 2 are defined in equation 
(1-2). 

In Section 4, we calculate diagram lc, and find that (in the lepton flavour basis 
where the sneutrino has no vev) it is approximately 

17113 ~ ~64^ mx - ( } 

Note that we neglect the supersymmetric partner of these diagrams, figure 3, 
because it is a small correction to the tree-level mass. As explained above, the GH 
diagram gives a mass to a different combination of neutrinos if Bj is not exactly 
aligned with \i j, that is Bj ^ Bjij. So in the basis where the sneutrino has no vev, 
we are only interested in the component of Bi which is orthogonal to Writing 



B i - ^r^fS = B\ = B l - B\, (3.10) 



we find 

ma m x . (3.11) 

The component of Bi parallel to /i; contributes to the mass of the neutrino that is 
massive at tree level. We also neglect the diagram ~ S^S^hf. This comes from mixing 
the neutrino with the up-type higgsino, h u , which has a Majorana mass from a loop 
like figure lb with tops and stops instead of bottoms/sbottoms. 

There is an additional contribution to the neutrino mass matrix in the case 
that H8b ~ S\. This comes from the diagram of figure la, with an L violating mass 
insertion at point VI. One unit of lepton number violation comes from a trilinear, and 
another from Bj, or m\ k . We use the minimisation condition for the Higgs/slepton 
potential: — tan (3 Bk = fn\ k (in the < v >= basis) [12], to get the neutrino mass 
estimate in table 1. 

3.3 Contributions when 5^ ^ 0, 5 B ^ 

In this section, we include L violating mass insertions on all possible propagators for 
diagrams with charged particles in the loops. Consider first the diagram of figure la. 
On the scalar propagator there can be a mass insertion 5b discussed at the end of 
the previous section, or 8^ (the latter appears in the off-diagonal mass term mixing 
the E c with L or H). Both external neutrino legs and the internal lepton line can 
have mass insertions 5^. Each of these possibilities is represented by a blob on the 
diagrams of figure I. The diagram is shorthand for a series of loops. The possibilities 
are listed in table 1, along with the corresponding estimates for m v . 

The diagrams must have two units of same-sign lepton number violation to gen- 
erate a neutrino mass. Diagrams with AL = +1 and AL = —1 would contribute to 
wavefunction renormalisation. In diagram la, this implies that we need one unit of 
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lepton number violation from I — III and one from IV -VIII. Schematically we expect 
terms of order 

(5„ + Sx + <V) (<V + 5 B + 5 X + <y (3.12) 
where the first parenthese corresponds to the blobs I III and the second to the blobs 
IV- VIII. We neglect for convenience in the table the m\ k insertions because these can 
be rewritten in terms of Bi using the minimisation conditions [12]. The contributions 
of order Sf, 5y, 5\5b have already been discussed, the corresponding estimates are 
given on the first, second and fifth lines of table 1, respectively. We briefly discuss 
here how we calculated the results presented in table 1 for the new diagrams which 
we consider. 

We first include the 5^ perturbations on the external lines at points I and VIII. 
These will be present for both A (figure la) and A' (figure lb) and correspond to 
expressions on lines 4 and 6 of table 1. We neglect the diagram with mass insertions 
on both legs because this will be proportional to the tree level mass ~ //j//j/m x . We 
can estimate the diagram with lepton number violation at points I and VII as follows. 
In the < Uk >= Vk — basis (k : 1..3), the incident neutrino Vi will be in the direction 
(vi) J = %A MJi /K- The mass insertion on the external leg is /Xj = v- X 1 ■ ft/m? cx S l , 
and allows the incident neutrino to turn into a Higgsino h u . In an arbitrary basis, 
where v j ^ is allowed, the mass insertion will be the /ijt>4 — Vi/j,4 oc \x ■ \ % ■ v. 
Continuing along the incident line, in the basis where the sneutrinos do not have 
vevs, the h u can turn into an hd via /i 4 , and at the vertex h d interacts with E c and £ 
via a Yukawa coupling. In an arbitrary basis, the neutralino after the mass insertion 
of fa or Vi will be an h u or a gaugino. These can respectively be turned into the linear 
combination of lj corresponding to the hd by /// or vj. These vectors are misaligned 
by an amount 5^, which we can neglect here because we have one unit of L violation 
from the mass insertion. We assume \v\ ~ \ ~ msusY- In column four of table 1, 
we therefore write the vertex as fii\' Ipq [fii\ IJk ] for the squark [charged] loop. The 
remainder of the diagram is the same as lines 1 and 2 of the table. 

Additional terms of order 8^(5^ + 5\) are obtained when we include the misalign- 
ment between jl and v at point V of diagram la, giving the contributions on lines 
7-9 of table I.. We can express 

p,= !—z-v + - 5 — v-X, 3.13 

v 2 mj 

or alternatively the part of (vrA + fj jR v u )\ RTk that is misaligned with respect to 
v R \ RTk can be written as 

pL R X Rpe v P v Q X^ se X STk 



(3.14) 



Finally there are corrections proportional to 5^ from /ij mass insertions on the 
internal lines of figure la. We assume vj\ JLl \ LMk VM = 5 tk mim k (because v — fi 
misalignment is higher order effect) to get the results listed in lines 9-13 of table 1. 
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If 5^ 7^ 0, the w mixes with the charged leptons, so there are loop diagrams with 
a gauge coupling at one end and a Yukawa coupling at the other, e.g. see figure Id. 
These gauge x yukawa diagrams can only appear if the gaugino mixes with the lepton 
on the internal line, so they are proportional to 5^. The other unit of L violation 
must be due to 5 M (at points I or VIII) or 5b (at point IV). These contributions are 
estimated in lines 14-16 of table 1. There is no contribution with a trilinear A at 
the vertex VII because it gives [m l J] ~ A lJ ' fc , and X ljk is anti-symmetric on indices ij 
whereas the mass matrix is symmetric. 



4. The GH loop in the MSSM mass eigenstate basis 

In this section, we calculate the one loop neutrino mass induced by slepton-higgs 
mixing (figure lc). We first do this using mass insertions the MSSM mass eigenstate 
basis, then compare to our previous exact (and basis independent) result in one 
generation [12]. The two calculations agree, but differ slightly from the original 
calculation of [7]. 

In the < v >= basis the two units of lepton number violation come from Bi and 
[ m L\1i (the crosses on the scalar propagators in figure lc), which mix the sneutrino 
with the Higgses. The scalars that leave the vertex where the neutrinos go in are 
sneutrinos. We can write [wii]^ = — tan f3Bi, using equation (3.8) in the basis where 
< u >= 0. 

The amplitude for this diagram connecting a neutrino ui to Uj is 
Aij = B i B~Y^ m Xc( Z a2 ~ Z al g'/g) 2 x 

Xa 

| (cos 2 a + 2 tan (3 cos a sin a + tan 2 (3 sin 2 a)/(m 2 , m 2 , rn 2 Xa ) 

+ (sin 2 a — 2 tan (3 cos a sin a + tan 2 (3 cos 2 a)I(m 2 H , m~, m 2 Xa ) 

- (cos 2 13 + 2 tan (3 cos (3 sin (3 + tan 2 (3 sin 2 (3)1 (m 2 A , m 2 , m\ a ) } . (4. 1) 

Recall {Xa} for a : 7. .4 are the MSSM neutralinos, and h, H, and A are the usual 
MSSM fields: 

h = cos aH* - sin aH% (4.2) 
H = sin aH^ + cos aHf (4.3) 
A = cos (3Hl - sin (4.4) 

(H^ and are the real and imaginary parts of the neutral components of the MSSM 
Higgs fields) and 



■ lgen 



(m 2 s ,ml,m 2 x J 



d A k 1 1 1 

(27r) 4 k 2 + m 2 k 2 + m 2 k 2 + m 2 k 2 + m 2 



Xa 
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In 



Xa 



{m? Xa - ml) (ml - m 2 ) ml 
mi , m 2 



(m 2 — m 2 )(m 2 —m 2 ) m 2 



ln^ . (4.5) 



This gives 



Q l 

m ^ = 4cos 2 /3 E( Z «2 ~ Z ^9'/9) 2 m Xa {l(m 2 h , m 2 , m 2 x J cos 2 (a - /3) 

+I(m 2 H ,ml,m 2 x J sin 2 {a- (3) - I(m 2 A ,ml,m 2 x j} , (4.6) 

which does not quite match Grossman and Haber's result, but goes to the same limits 
when various masses become large. We can check equation (4.6) by expanding the 
exact result of [12]. The result of [12] is for one lepton generation only; we assume 
here that the sneutrinos are degenerate, so we can compare to this calculation. From 
[12] we have: 

= £2 E m Xa(Z a2 - Z al g'/g) 2 J2(v • s n ) 2 e n B (0, M 2 : m 2 x J (4.7) 

D47r a:7.A n 

where e n is +1 for s n a CP-even scalar, and -1 for n a CP-odd. Note that in an ]fl p 
theory, the sneutrinos split into a CP-even and a CP-odd scalar, that are not mass 
degenerate, v is the lepton direction corresponding to the incident neutrino, and s n 
is the direction in Higgs-slepton space corresponding to the mass eigenstate s n . B 
is a Passarino-Veltman function: 

£? (0, M 2 ,m 2 ) = -167r 2 *lim ' <M ' 1 




(4.8) 

There are divergent and scale-dependent contributions to B in addition to I(M 2 , m x ); 
however these cancel in the sum over scalars Sj in equation (4.7). 

To compare this result to equation (4.6), we need to expand the mixing angles 
(v-s n ) and the masses M 2 in £>j. We define the neutrino direction to be z>j = v-X jm^ 
(this is the definition of the charged lepton direction we have been using all through 
the paper). 

As Ifip^ 0, we make the following identifications: 

s 5 = h -> h 
S4 = hi 

s 3 = h 3 -+ s 8 = A 3 -> v\ (4.9) 
s 2 = h 2 

Si = /li 



H 


S9 


= A,- 






S8 


= A 3 - 




v R 


S7 


= A 2 - 




~ V R 


«6 


= M- 





13 



where h n is CP-even, A n is CP-odd, v R is the real part of the sneutrino, and the 
right hand side of the arrow are MSSM fields. 

The diagram that generates a Majorana mass for a neutrino must contain two 
units of lepton number violation. There is a contribution from the mixing angles 
between the neutrinos and the s n ,n : 9..1, and from the hk — A^ mass differences for 
k : 3..1. We consider first the contributions from /is,/i4 and At. The mixing angles 
(z>j • s n ) for s n = h 5l h^ or A 4 will be proportional to B i: so I/{ p corrections to the 
/i 5 , /i 4 or A 4 masses can be neglected in expanding equation (4.7). They would be a 
higher order effect. We find, for instance, 

' £>i(cos a + sin a tan 3) , 

Vi-h b = 2 ' ( 4 - 10 ) 

ml - ml. 

which substituted into (4.7) gives the contribution of the last line of equation (4.5) 
to the complete result. 

The contribution from /i 3 and A 3 ( scalars that are mostly u T ) can be written 



As E m Xa (Z a2 - Z al g'/gf {/(M 2 3 ,m 2 J - 

a:7.A 

-[(y T ■ kf + (u T ■ Kf - (u T ■ A^IiM^ml)} . (4.11) 

The first line is the contribution from the v R — v 1 mass difference, and the second line 
is from the mixing angle: {i' T -h 3 ) 2 =l — (£> T -hc > ) 2 —{0 T -h i ) 2 . The third term of equation 
(4.11) gives contribution from the second last term of equation (4.5) to the complete 
result, and the first two terms of equation (4.11) give the contribution from the first 
two terms of equation (4.5). This can be seen by writing m 2 h . A = m? T + Am^/2, and 

^2 ( B 2 (cosa + sin a tan/?) 2 B 2 (sin a — cos a tan 8) 2 



Ut \ m\ — ml m\ — ml 

B 2 (cos3 — sin 3 tand) 2 \ 
+ — —2 ~2 — • 4 - 12 

m A ~ m v J 

The result in [7] corresponds to this contribution. 

We can make a similar expansion for ■ h 2 and z> M • A 2 , and for v e assuming as 
here that there is no flavour violation amoung the sneutrinos, e.g. [y T ■ A 2 ) = 0, etc. 



5. Discussion 

We now return to our introductory suggestion that neglecting IfL p masses in the 
loops creates confusion about basis. Consider the superpotential for a one-generation 
model. This could have two ]f{ p parameters: A' or fi 3 . As discussed after equation 
(1.3), one can define the Higgs such that /x 3 = or such that A' = 0. In the basis 
where 113 = 0, there is the usual loop neutrino mass corresponding to diagram lb 
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with ]fi, p dX the vertices II and VII. Now consider the basis where A' = 0. If the mass 
/i3 is neglected in the loops, then it appears that there is no loop neutrino mass. 
This is perplexing, because we know that the tree-level mass is proportional to the 
misalignment between the Higgs-slepton vev and (/i 4 , If we are in a model where 
jl || v, the basis rotation cannot have moved the neutrino mass from loop to tree level. 
Thus, it is unclear where did the contribution to the mass go. If we include the loop 
lb with ]fi p at points I and VIII, this confusion is removed. 

We conclude that the basis in which neutrino mases are calculated is irrelevant, 
providing that all the contributions are included. We believe that "basis-related" 
confusion comes from neglecting the ]ft p masses in computing loop neutrino masses: 
the magnitude of the contribution being neglected depends on the basis choice. 

We therefore sit in the basis where the sneutrino has no vev, and consider the 
size of different contributions and the structure of the neutrino mass matrix in lepton 
flavour space. The tree-level mass contribution to the neutrino mass matrix is 

Kk ~ (5.1) 

m x 

in all cases. 

In case A, where we neglect ^masses in the loops, the loop contributions from 
the squark/quark and lepton/slepton loops are respectively of order 

Note that the quark basis can be chosen to simultaneously diagonalise the A' and 
down-type Yukawa couplings, but that the lepton basis that diagonalises the lepton 
mass matrix X IJk vi does not neccessarily diagonalise X l£k on indices £k. 

In case B, we include the soft masses {Bi} and {m^} in the GH diagram (figure 
lc) which contributes 

where and Bj are defined in equation (3.10), and we have neglected the term of 
order BfBj because it is aligned with the tree-level mass of equation (5.1). 

In case C, we include all the lfi p masses in perturbation theory. All the terms 
listed in column 5 of table 1 will contribute. Those from diagram Id, which involves 
a gauge coupling, are potentially largest: 

M y ~i6^K^ 2 + < m » 2 ]# (5 ' 4) 

Note that this is not aligned with the tree-level mass due to the inclusion of lepton 
masses. There are additional contributions of order h A jiijj,j/\ji\ 2 or h A Bijj,j/ \\n\B) 
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from mass insertions on diagrams la and lb. (h is some Yukawa coupling.) These 
are listed in table 1. 

To summarize we have estimated one-loop contributions from R-parity violating 
parameters to the neutrino mass matrix in the context of the MSSM. We included 
the fi p bilinears in the mass insertion approximation, which introduces new loops 
and new contributions to the canonical graphs considered in the literature. For the 
neutral graph arising from slepton-Higgs mixing, our mass insertion results are in 
agreement with the perturbative expansion of the exact result. 
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Figure 1: Schematic representation of one- loop diagrams contributing to neutrino masses, 
in a Lagrangian basis. The blobs indicate possible positions for Ifip interactions, which can 
be trilinears (at positions II and VII) or mass insertions. The misalignment between jl 
and v allows a mass insertion on the lepton/higgsino lines (at points I, III, or VIII) and at 
the ^4-term on the scalar line (position V). The soft Ijl v masses appear as mass insertions 
at positions VI and IV on the scalar line. Figure a) is the charged loop with trilinear 
couplings A (or h e ) at the vertices. Figure b) is the coloured loop with trilinear A' or 
yukawa couplings. Figure c) is the neutral loop with two gauge couplings (this diagram 
is drawn in the MSSM mass eigenstate basis.), and figure d) is the charged loop with one 
gauge and and a Yukawa coupling. This diagram is relevant if gauginos mix with charged 
leptons — that is if 5^ ^ 0. 
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z 



Figure 2: The Grossman-Haber contribution in interaction eigenstate basis. The scalars 
leaving the neutrino vertices are sneutrinos. The crosses on the scalar line are Bi, Bj, m|j 
or m\y 




Figure 3: Feynman diagram showing the gauge loop correction to the tree-level neutrino 
mass in MSSM neutralino mass eigenstate basis. This is the supersymmetric partner of 
the GH loop. The small "x"s on the internal fermion line correspond to the interaction 
\iiVih u . The large X is a neutralino mass. There would be additional diagrams in a basis 
where < v >^ 0. 
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Table 1: Estimated contributions to \m v ] 13 from all the diagrams. In the 
second two columns is the label of the diagram of figure 1, and the position 
on the diagram of the two AL = 1 interactions. Column four is the "basis 
independent" combination of coupling constants for the diagram (which must 
be symmetrised on i H j), and column five is the estimated contribution 
to the neutrino mass matrix (in the "flavour basis" where < i> >— and 
i)i = v ■ X z /mf). All indices other than i and j are summed. 



